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or the past several years, 
the creation of computer 
generated fractal images 
hhas become a popular pas- 
time with many program- 
mers. In particular, the Mandelbrot Set is 
probably one of the most studied fractals 
today. 


What Is A Fractal? 


In simple terms, a fractal is a compli- 
cated shape which continues to remain 
complicated no matter how much you 
magnify it, For example, consider the 
famous Snowflake Curve, whichis gener- 
ated as follows: start with an equilateral 
triangle. Trisect each side of the triangle, 
and on each of the middle thirds erect an 
equilateral triangle pointing outward. 
Erase the parts common to the new and 
the old triangles. The result looks like the 
outline of a "Star of David.” Next, trisect 
each of these twelve lines and again upon 
each middle third erect an equilateral tri- 
angle pointing outward. Continue this 
process forever. The result is a fractal 
image resembling a snowflake. 

Another interesting fractal is the 
Dragon Curve. To generate a Dragon 
Curve, simply select two points and con- 
neetthem witha jagged path consisting of 
straight line segments. This path then 
becomes the “rule by which you move 
from any point to any other point. In 
particular, the jagged path you specified 
is itself a series of points connected by 
line segments. Therefore, each of these 
line segments can be replaced with a 
smaller copy of the original jagged path. 
Continuing this process recursively, 
yields a Dragon Curve. 

Fractals need not be limited to planar 
images. For example, consider the fol- 
owing process: start with a large wooden 
cube with faces facing north, south, east, 
‘west, up, and down. Now, saw this cube 
into a 2x2x2 set of smaller “cubies” (by 
making three perpendicular cuts). Next, 
throw away the cubies in the top north- 
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‘west corner, the top southeast comer, the 
bottom northeast comer, and the bottom 
southwest comer. What remains are four 
cubies joined together along their edges. 

Repeat this process for each of these four 
remaining cubies. The result will be six- 
teen “sub-cubies” joined together along 
their edges. Next, saw each of these six- 
teen sub-cubies... The result is a fractal 
structure that can probably best be de- 
scribed asa “filigreed tetrahedron” stand- 
ing on one of its edges. 

This tetrahedron structure has the in- 
teresting property of looking like a solid 
square when viewed from each of the six 
directions (north, east, etc.) Furthermore, 
every sub-cubie is visible from the six 
directions (ie., no two sub-cubies ever lie 
‘on the same line of sight). Therefore, you 
could take six different images and cut 
them up into little squares, and paste them 
on the faces of the cubies in such a way 
that a different one of the six original 
images is completely visible when you 
view the structure straight-on from each 








of the six directions. (actually made one 
of these things out of 256 wooden build- 
ing blocks.) 


The Mandelbrot Set 


‘The American public first became 
aware of the Mandelbrot Set when it 
appeared as the cover story of the August 
1985 issue of Scientific American. In his 
‘monthly column titled “Computer Rec- 
reations,” A. K. Dewdney described how 
to generate this “most complex object in 
mathematics.” 

In the complex plane, the Mandelbrot. 
Sets defined to be the set of points, c, for 
which the process: 

24:5 2,'+62,=0 (Equation #1) 
remains bounded forall values of n. (This 
is represented by the black region in Fi 
ure 1.) However, when we generate 
Mandelbrot images, we are generally not 
as interested in the Mandelbrot Set itself 
as we are in the points that lie just outside 
of the set (ie., those values of ¢ for which, 
equation #1 is almost bounded, but which 
starts to diverge for sufficiently large n). 
If we let: 


z=a+biand 








=xeyl 


where i is equal to the square root of -1, 
then equation #1 can be expressed as: 


242 82- DP 4x 
(Equations #2) 


b,=2a,b,+¥ 
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Generating Mandelbrot 
Images 


Using equations #2, we can write the 
following little BASIC subroutine: 


1000 
1010 
4020 
1030 


4040 


1050 


This subroutine takes an X and Y 


TAT SY:A=T 
41 :IFN> 255 THEN 





IFAN24B42>4 THEN 


coordinate (the real and imaginary com- 
ponents of c), and returns the valueN for 
which equation #1 start to diverge. By 
using a pair of nested FOR/NEXT loops 
to“raster scan” the values of X and Y over 
aregionof interest, and using the returned 
value of N to select the pixel color, a 
Mandelbrot image can be generated. 
However, I do not recommend using this 
subroutine for generating images. Even 
‘on an Apple 11GS running in fast mode, 
this litle BASIC subroutine can take four 
orfive hours to generate even a40x40.Lo- 
Res image! 

‘A faster way to generate Mandelbrot 
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pixels is to use assembly language. Such 
‘a subroutine is shown in Listing 1. (in a 
benchmark test, this assembly language 
subroutine ran about 30 times faster than 
the previous BASIC subroutine.) 

Subroutine MANDEL uses a special 
purpose multiplying subroutine (MULT) 
shown in Listing 2, In addition, it refer- 
ences the familiar SAVEO and UN- 
‘SAVEO entries that! presented backin the 
January issue of Call -A.P.P-L.E. Since 
it’s been more than six months since they 
have appeared, I am republishing them 
this month (Listing 3). 

‘Subroutine MANDEL has four call- 
ing parameters. MANDEL-12 specifies 
the timeout value for the pixel (i.e., the 
number of cycles you want to spend in the 
loop). If equation 1 has not diverged by 
this time, the pixel is assumed to be in the 
Mandelbrot Set. (The default value for 
this parameter is $0100=256.) The two 
bytes at MANDEL-10 contain the pixel 
value returned by the subroutine. This 
pixel value is related to the n in equation 
#1, except that itruns backwards from the 
timeout value. In other words, a returned 
value of zero means the pixel lies in the 
Mandelbrot Set, while a large value 
theans that the pixel lies far away from the 
Mandelbrot Set. 

‘The real and imaginary components 
(ie., X coordinate and Y coordinate) of 
the pixel are passed in MANDEL-8 and 
MANDEL-4, respectively. Each of these 
coordinates is a 4-byte fixed point num- 
ber. The integer portion scontained in the 
most significant byte while the fraction 
portion is contained in the least signifi- 
cant three bytes. Negative numbers are 
expressed by complimenting all the bits in 
the four bytes. For example, the: value 
25 would be expressed as $01400000, 
while the value-1.25 would be expressed 
as SFEBFFFFF. 

Notice that subroutine MULT uses a 
dozen or so locations on page zero with- 
out first saving those Toeations. Is this a 
violation of our Rule 2? (See “Let's Start 
Publishing Usable Assembly Language 
Code,” Call-A.P-P LE. December 1986.) 
If subroutine MULT had been intended 
for general purpose use, then it would 
indeed be violating Rule 2. But MULT is 
a special subroutine that is intended to be 
called only from subroutine MANDEL. 
(Itis what the PACK documentation re- 
fers to as a “cooperative” subroutine.) In 
fact, subroutine MANDEL first invokes 
subroutine SAVEO before subroutine 
MULT ever gets called. We can think of 
subroutine MULT as merely being an 
extension of subroutine MANDEL which 
was broken off into a separate source 
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listing in order to keep the source listing 
for MANDEL from becoming too long. 


Generating Lo-Res 
Mandelbrot Images 

Listing 4 is an Applesoft program for 
generating Lo-Res Mandelbrot images. 
(The program assumes that the MAN- 
DEL module has been loaded at memory 
location $8000, so that the MANDEL, 
entry is located at $800C.) 

‘The program asks for three input para- 
meters: the X coordinate, the Y coordi- 
nate, and the magnification. The X and Y 
‘coordinates are merely the real and imagi- 
nary components of the region you wish 
to view. (The image will be centered on 
these coordinates.) The magnification 
specifies, as you would suspect, how big 
you want the image to be. To test out the 
program, try the following: 

X COORDINATE: -0.7456185 

Y COORDINATE: +9.1050935 

MAGNIFICATION: 600 


The resulting image is shown in Fig- 
ure 2. It took about seven and a half 
minutes to generate this image on an 
Apple Ils. It should therefore take about 
15 or 20 minutes on a regular Apple IL 

‘When the program finishes, it beeps the 
speaker every few seconds to let you 
know it’s done. 


Generating Hi-Res 
Mandelbrot Images 


Listing 5 showsan Applesoft program 
for displaying Mandelbrot images on the 
Hi-Res sereen. (Again, the program as- 
sumes that the MANDEL subroutine 
module has been loaded at $8000.) Notice 
that lines 10-1080 of the listing are ex- 
actly identical to the beginning of the 
program shown in Listing 4. Therefore, 
you can simply delete the last few lines 
from Listing 4 and replace them with the 
new lines (2000-3040) of Listing 5. 

Totest this program, try specifying the 
following: 

X COORDINATE: -0.66666666 

Y COORDINATE: 0 

MAGNIFICATION: 0.4 





(es, you can use magnifications less 
than one.) The result is the image of the 
entire Mandelbrot Set that we saw in 
Figure 1. 


Limitations 


Those readers who are intimately 
familiar with the Mandelbrot Set may 
have noticed that the contour lines (equi- 
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potential lines) in Figure 1 are not quite 
correctas you getfurtherand further away. 
from the Mandelbrot Set. This is because 
our criterion for divergence was to test for 
the size of z exceeding the value 2. This 
was the value that A. K. Dewdney sug- 
gested in his Scientific American article. 
To examine points far away from the 
Mandelbrot Set we need to use a value 
much larger than 2. But as long as we stay 
close to the set (the region that’s really the 
‘most interesting anyway), the value of 2 
‘works just fine. 

‘There is also a limit to how much 
‘magnification we can use when generat- 








ing an image. This is due to the finite 
precision with which the calculations are 
being performed. Therefore magnifica- 
tion should not exceed approximately 
25000. At very high magnifications, 
round off errors start to become a serious 
problem. For example, Figure 3. was 
‘generated with a magnification of 15000. 
Ifyou look closely, you can see that some 
of the pixels are already beginning to 
show signs of error. 





Creating “Pretty Pictures” 
“Mandelbrot images representan nte- 





FIGURE 4 
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esting blend of technology and art. Even 
though Mandelbrot data can be generated 
almost mechanically, the displaying of 
that data as an image is far from being cut 
and dried. Some of you may have already 
tried to generate Mandelbrot images in 
the past but could never quite get them to 
look like the pictures that are published in 
books and magazines. 

‘Once you have generated the Mandel- 
‘rot data, you must decide how you want 
to display it. For example, Figure 4 was 
generated by the program shown in List- 
ing 5 with the following parameters: 





X COORDINATE: 40.358 
Y COORDINATE: 40.643 
MAGNIFICATION: 35 


The image shown in Figure 5 was 
created from exactly the same data, ex- 
cept that this data was processed a little 
differently. Instead of simply assigning 
each retumed pixel value a color (as was 
done in Figure 4), the data was high-pass 
filtered. In other words, the raw value for 
‘each non-zero pixel was compared with 
the raw value of the pixel to its immediate 
left and with the raw value of the pixel 
immediately below it. If either of these 
two neighboring pixels differed from the 
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current pixel by no more than plus or 
minus 1, the current pixel was assigned 
the color black. This eliminates all the 
contour lines and displays only those 
pixels that are close to the Mandelbrot Set 
(ie., those pixels that come from struc- 
tures that are too small to be seen at the 
current magnification) 

A similar but slightly different tech- 
nique was used to create the super Hi-Res 
image shown in Figure 6. Here, the con- 
tour lines were allowed to remain, Butthe 
“noisy” pixels (the ones that are not re- 
lated totheirneighbors) were assigned the 
color white. 

The super Hi-Res images shown in 
Figures 7 and 8 were both created by yet, 
another technique. (Figure 7 appears on 
the cover.) There, the color assigned to 
each pixel was determined by its coarse 
value rather than by its fine value. In other 
words, any pixel whose value was greater 
than zero, but less than a certain level was 
assigned one color. Those pixels that 
‘were greater than that level, but less than 
a higher level were assigned a slightly 
different color, and so on. By judiciously 
choosing shades of yellow and orange to 
match the flame-like structure of the 
image, we obtained the fiery picture 
shown in Figure 7. Similarly, by assign- 
ing shades of blue, white, and purple to 
the pixels from yet a different region, we 
obtained the “electrified” image of the 
very tiny Mandelbrot satellite shown in 
Figure 8. (The coordinates of this image 
are: -0.1634, +1.0261, with a magnifica- 
tion of 3500. In addition, the timeout 
value was increased to 333.) 











Things To Do On Your Own 


In part IV of this series (see Call - 
APPLE,, April 1987), I published sub- 
routineRWBLOCK, amodule foraccess- 
ing the 64K auxiliary memory. AS one 
application of that subroutine I presented 
a simple frame Hi-Res movie. You 
might like to modify the program shown 
in Listing 4 touse subroutine RWBLOCK 
and automatically save a sequence of Lo- 
Res images of the Mandelbrot Set, and 
then play them back in the form of a Lo- 
Res movie. However, be prepared to lock 
your computer away for the weekend 
since such a movie will take several days 
to generate! ([ have made several such 
“Fractal Flickers.” By generating a se- 
quence of Lo-Res images, each with a 
slightly greater magnification than the 
previous one, I was actually able to zoom. 
down into different region of the Mandel- 
brot Set.) 

‘The subroutine library in volume I of 
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the PACK contains subroutines that let 
you plot pixels on the double Hi-Res 
‘graphics screen. Those readers who have 
access to the PACK might want to modify 
the program in Listing 5 to generate 
double Hi-Res Mandelbrot images. (For 
more information on the PACK, contact: 
Interactive Arts, 2715 Porter Steet, 
Soquel CA 95073, telephone: (408) 475- 
7047. 

‘Subroutines MANDEL and MULT 
work with 4-byte coordinate values. 
Consequently, the maximum  magni- 
fication that can be achieved is ap- 
proximately 25000. You might want to 
try to increase the magnification range by 
1g Listings 1 and 2to incorporate 
more bytes of precision in their calcula- 
tions. In principle, there is no limit to the 
resolution that can be achieved. By incor- 
porating enough bytes of precision, you 
can generate images that let you “boldly 
‘go where no man has gone before.” 

‘The Mandelbrot Set is not the only set 
‘of images that can be generated by the 
equation: 

















If, instead of starting with: 
220 


and raster scanning c, we let ¢ remain 
constant and raster scan different initial 
value for z, we obtain images of the so- 
called “Julia Sets.” While the Julia Sets 
are related to the Mandelbrot Set, they are 
somewhat different and interesting in 
their own right. 


Conclusion 


Generating Mandelbrot images is a 
fun pastime. (Besides, it gives your com- 
puter something to do while you're 
asleep!) But, even using machine lan- 
‘guage, it can take several hours to gener- 
ate an entire image. Therefore, some en- 
terprising software developer might 
someday come out with a Mandelbrot 
“coloring book” program consisting of a 
disk full of raw Mandelbrot data for sev- 
eral different regions, along with a set of 
user-friendly image processing programs 
for displaying that data with different 
filters, colors, etc. Such a package would 
be particularly fun to play with on the 
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mode. 

Even though each pixel of the Mandel- 
brot Set is generated independently, the 
pixels are still calculated one at a time. If 
you could have several computers. all 
‘working together on different parts of the 
same image, the entire image could be 
created just that much faster. With this as 
a motivation, several of my colleagues 
and I are currently designing an inexpen- 
sive, general purpose “super-computer” 
consisting of up to 256 6502 micropro- 
cessors all working in parallel and run- 
ning at 3Mhz. The interface for this de- 
vice will plug into one of the Apple's 
expansion slots so that it will be accessed. 
as a peripheral. When programmed to 
‘generate Mandelbrot data, the Apple I 
will then simply become a graphics termi- 
nal for displaying the resulting data. Allof 
the “number crunching” will be done 
inside of the parallel processor machine. 
When fully functional, we expect to be 
able to generate complete Mandelbrot 
images in less than a minute! 


Please rate this articie on the Reader 
service Card by circling: 
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ES TAX 
0 Dex 
41 PULER PLA 
2 STA $00X 
ct Dex 
108 a BL PULLER + 
4S UNSAVED PLA 
6 Tay 
a PLA 
4s TAK 
“9 RTS 
aos 
003. St TEMP. Og ag 
So 
LISTING 4 
40 INPUT X COORDINATE: “REAL 
20. INPUT COORDINATE: “IMAG 
90. INPUT “MAGNIFICATION: “MAG 
40 SZE-1/MAG 
50 SUB=87 4006+ 12 
90 GOsuB 2000 
400 XO AEAL- SIZE/2+ SIZE /(2*NX) 
150 OX=SIZE/NK 
200 DY =DX/ ASPECT. 
250. YY = IMAG -DY “NY /2+DY/2 
900 FORY=1 TONY 


sto 


soeee8 


g2es88 


Vir YY; GOSUB 1000 
POKE SUB- 4,8: POKE SUB -3B 

POKE SUB. 2/6: POKE SUB- 1.0 

X= x0 

FORX-=1 TON 

= 0%: GOSUB 1000 

POKE SUB -8.A: POKE SUB-7,.8 

POKE SUB -6.C: POKE SUB- 5,0 

‘CALL SUB 

2 PEEK (SUB- 10) +256" PEEK (SUB -9) 
‘Gosue s000 

2X = XX + DX: NEXT X 

YW YY + DY: NEXTY 


1010 


2H0 
2120 
2190 


2010 


S88eseeseus 


e8eee82 


Sugasss 


(CALL 64477: REM BELL 
FOR Z~ 1 TO 10000: NEXTZ 
‘GOTO 600 

$= SGN().V~ ABS (V) 

Da INTV=V-D 
Vn256°V.C= INT (WIV =V-C 
Ve 286°VB = INT(VIV=V-B 
Vn286°VAm INT(VV AVA 
IFS>OTHEN RETURN 
0=255-D:c-255-C 

Ba 255-BA~255-A 
RETURN 

Nx= 4ONY = 48 
ASPECT = 153 
TOL(i) = -TELY) «9 
TBL(2) = 13:7BL(a) = 12 
BLS) = 6TBL6)=3 

POKE SUB 120 

POKE SUB 11,1 

GR: POKE - 169020 
FORY ~ 40.70 47 
HUNO39ATY 
‘NEXT ¥: RETURN 

IFZ=0 THEN RETURN 
ZaZ-6" NTZ/8) 


SESPRTREY: Herurn 











LISTING 5 


ae 
SUT + S3RBNATE: Sums 
INPUT MAGNIFICATION: AG 
SIZE» 1/MAG 

SUB = 8 408+ 12 

osv 2000 

XO REAL-SZE/2.+SZE7@*N) 
Dx SE/NK 

DY=OX/ASPECT 
YWaIMAG- DY "NY/2«D¥/2 
FORY-1 TONY 

V=Y¥: GOSue 1000 





Ya 4 DY: NEXTY 
‘CALL 64477: REM BELL 
‘ORZm 1 TO 10000: NEXT Z 
goTo.800 

‘Sm SGN (VV = ABS (¥) 

Do INT(WV~V-D 

Vn 258°V.C= INTVV=V-C 
V=256° VB = INT (VV <V-B. 
Vn256°VA= INT (VV =V-A 
18> OTHEN RETURN 
Dn255-D:0=255-C 

82255 -BA=255-A 
RETURN 

NX H40NY » 192 
ASPECT = 1.75 

POKE SUB- 120 

POKE SUB-1.4 

HGR: RETURN 

IFZ =O THEN RETURN 
Za2-3" INTZ/9) 
HOOLOHRZ +5 

HPLOT2* 0c-1),182-¥ 
HPLOT2°X- 1,192- ¥:RETURN 


Call-AP.PLE® August 1967 


